Gate tunable quantum transport in double layer graphene 
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We analyze the effect of screening provided by the additional graphene layer in double layer 
graphene heterostructures (DLGs) on transport characteristics of DLG devices in the metallic 
regime. The effect of gate-tunable charge density in the additional layer is two-fold: it provides 
screening of the long-range potential of charged defects in the system, and screens out Coulomb in- 
teractions between charge carriers. We find that the efficiency of defect charge screening is strongly 
dependent on the concentration and location of defects within the DLG. In particular, only a moder- 
ate suppression of electron-hole puddles around the Dirac point induced by the high concentration 
of remote impurities in the silicon oxide substrate could be achieved. A stronger effect is found 
on the elastic relaxation rate due to charged defects resulting in mobility strongly dependent on 
the electron denisty in the additional layer of DLG. We find that the quantum interference cor- 
rection to the resistivity of graphene is also strongly affected by screening in DLG. In particular, 
the dephasing rate is strongly suppressed by the additional screening that supresses the amplitude 
of electron-electron interaction and reduces the diffusion time that electrons spend in proximity of 
each other. The latter effect combined with screening of elastic relaxation rates results in a peculiar 
gate tunable weak-localization magnetoresistance and quantum correction to resistivity. We propose 
suitable experiments to test our theory and discuss the possible relevance of our results to exisiting 
data. 



I. INTRODUCTION 

Gate tunable Anderson localization of Dirac electrons 
is inaccessible in standard field effect heterostructures 
due to the strong disorder induced inhomogeneity of the 
local doping level (electron-hole puddles )^~— ubiquitous in 
the materials that show Dirac-like low energy spectrum, 
such as graphene^ or topological insulator surfaces^. Ba- 
sically, the low-density localization behavior of graphene 
is inaccessible experimentally becasue the Dirac point it- 
self is inaccessible due to the formation of electron-hole 
puddles around the charge neutrality point. However, 
very recently, a novel double layer graphene (DLG) het- 
erostructure containing two graphene layers separated by 
an insulator allowed using additional screening effect of 
the second graphene layer to access a metal-insulator 
transition regime. Anomalously large resistnace, p 3> 
/i/(4e 2 ), in combination with insulating temperature de- 
pendence was observed in these experiments suggestive of 
a metal insulator transition (MIT)S. The physical nature 
of this behavior is currently a topic of a debate^^. In 
particluar, whether the transport data of Ref . represent 
an effective low-density high-temperature semi-classical 
resistivity^ or a low temperature strong Anderson local- 
ization crossover behavior— is unclear at this stage. 

On the one hand supressing inhomogeneity allows ac- 
cess to the low-density regime kp£ < 1 in which Ander- 
son physics is expected to dominate (here fc^ is the Fermi 
wavevector and I is the mean free path). In this regime, 
quantum interference of the two flavors (due to two val- 
leys) of chiral Dirac charge carriers in graphene may re- 
sult in the insulating behavior in the case of sufficient 
mixing of the flavors by atomic scale disorder—. On the 
other hand, the experiments of Ref. d were restricted to 
a fairly high temperature regime (10-100K) where quan- 



tum interference effects may be suppressed due to short 
phase breaking length. The observed MIT is also in con- 
trast with more recent transport measurements on ultra 
high quality suspended devices 1 ^ in which an extremely 
low density inhomogeneity is achieved, 5n ~ 10 8 cm~ 2 , 
nevertheless no MIT is detected. Moreover, an Anderson 
insulator would be characterized by the resistivity grow- 
ing exponentially with decreasing temperature whereas 
the data 9 demonstrates only a roughly power law growth. 
Insulating power law temperature dependence (resistiv- 
ity growing with decresing temperature) is not unsual 
near charge neutrality in graphene (see Fig. 2 in Ref. 
and measurements on other lo w 12 i 13 and high 11,14 mobil- 
ity samples). This behavior is explained by the combined 
effect of temperature dependent screening 1 ^ thermal ex- 
citation of electrons from the valence band and thermally 
activated hopping of electrons over the potential barriers 
between electron-hole puddles 1 ^. Nevertheless, very high 
values of resistence reported in Ref. d, p 3> h/(4e 2 ) sug- 
gest a novel behavior in this system possibly associated 
with MIT. However, an alternative explanation for the 
observed behavior was recently suggeste d 10 ! 16 , in which 
case the anomalous resistance is explained by a Boltz- 
mann transport effect combined with strong suppression 
of inhomogeneity in DL G 10 ' 16 . The latter work is phe- 
nomenological and relies on the assumptions of the qua- 
siclassical transport formalism and the phase breaking 
length being shorter than the elastic scattering length. 
Therefore a more detailed analysis of quantum effects 
is required to justify the applicability of the latter ap- 
proach. One way to make progress in the understanding 
of the observed behavior of the resistivity in these new 
devices is to analyze the metallic regime where perturba- 
tion expansion in disorder strength may be applied. Such 
theoretical analysis of the quantum interference effect in 
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the higher-density weak localization regime of the DLG 
system, which must be a precursor to any low-density 
strong localization crossover phenomenon, is one of the 
main goals of our work. 

In this paper we consider the effect of gate tunable 
screening provided by the additional layer in DLGs on 
transport characteristics of these devices. Screening af- 
fects both the potential of charged defects and Coulomb 
interactions of charge carriers. Taking into account both 
effects we analyze the classical and quantum parts of the 
resistivity of graphene in DLG geometry in the metal- 
lic regime kpt 3> 1. We consider only low temperature 
regime such that elastic scattering by the disorder limits 
the transport characteristics and all inealstic scattering 
effects may be included perturbatively. 

Below we outline the results of the following analysis: 
(i) We develop a framework describing the screening ef- 
fect of the two graphene layers in DLG on the Coulomb 
potential of charged defects in DLG hctcrostructure. (ii) 
We generalize the self-consistent theory of disorder in- 
duced electron density fluctuations at neutrality point of 
graphene^ (electron-hole puddles) and use it to estimate 
the screening effect of DLG on the amplitude of the elec- 
tron density fluctuations, (iii) We analyze the effect of 
screening in DLG on the elastic scattering rate due to 
charged defects that determines the extent of gate tun- 
ability of mobility in DLG. (iv) We analyze the effect 
of the screening in DLG on the dephasing rate due to 
inelastic electron-electron collisions, (v) We analyze the 
combined effect of screening on the weak localization cor- 
rection to the resistivity (vi) We predict a peculiar gate 
tunable magneto-resistance in the high quality DLG stru- 
tures with mixed elastic scattering mechanisms. 

The present work is limited to the perturbative metal- 
lic regime and cannot make any definitive conclusions 
about the regime of the low density studied in the ex- 
periments of Ref. [^. However, the theoretical frame- 
work developed below provides a basis for analysis of 
a more detailed gate tunable quantum transport mea- 
surements which would allow characterization of various 
microscopic scattering mechanisms that determine the 
behavior of the system in the low density regime. In 
particular, we believe that it is imperative that we un- 
derstand the high-density metallic regime of the DLG 
transport in some depth before trying to understand the 
low-density crossover MIT regime since a systematic per- 
turbative theory is available in the higher density regime 
whereas the low-density regime is inaccessible to analyt- 
ical theory. We hope that our work would motivate ex- 
perimental work focusing on the metallic regime of DLG 
in order to investigate conceptually important question 
of Anderson localization and metal-insulator transition 
in graphene. 

We consider a DLG consisting of two layers of 
graphene, "studied" and "control" , see Fig. [TJ which are 
used for resistivity measurements and to provide the ad- 
ditional screening effect, respectively. The two layers are 
separated by a thin layer, d > 4nm, of insulating hexago- 



nal boron nitride (hBN) that is thick enough to suppress 
direct tunneling between the layers. The whole strucure 
is separated by a thick layer, do ~ 20nm, of hBN from 
the standard Si02/Si substrate. The combination of the 
top and bottom gates in this heterostructure allows in- 
dependent "control" of the electron densities n c and n s 
in the "control" and "studied" layers respectively. This 
setup is motivated by the design of heterostructures in 
Ref. y. We assume that the mobility of the "studied" 
layer is limited by elastic scattering of charged impuri- 
ties in accordance with the data£. 

The additional screening effect of the "control" layer 
is expected to suppress the electron-hole puddles in the 
"studied" layer. We find that the effectiveness of the 
"control" layer screening strongly depends on the posi- 
tion of impurities in the DLG. It is efficient for the rel- 
atively small concentration of charged defects on top of 
the "control" layer, whereas it leads to only a moderate 
suppression of the inhomogeneity induced by charges in 
close proximity of the "studied" layer and by the remote 
charged defects in Si02 substrate. The latter inefficiency 
stems from the very high densities of charged defects in 
the Si02 substrate, nt > 10 12 cto -2 , that would require 
screening charge density in the "control" layer several 
times rii which is beyond the experimentally accessible 
range. Moreover, the nonlinear screening effects in the 
electron-hole puddle regime in the "studied" layer work 
against the screening effect of the "control" layer making 
the screening less efficient. The situation is slightly dif- 
ferent in the case of elastic relaxation rate due to charged 
defects. The screening by "control" layer is more efficient 
in this case and leads to the suppresion of the classical 
part of the resistivity near the Dirac point of the "stud- 
ied" layer. The resulting behavior of the resistivity as a 
function of the density in the "control" layer is a result 
of a competition of the suppression of electron-hole pud- 
dles that increases the resistivity and suppression of the 
elastic relaxation rate that decreases the resistivity. 

Below we also demonstrate that the ratio of momen- 
tum relaxation rate to the transport rate varies substan- 
tially with the distance between graphene and charged 
defects. Therefore measurement of the two elastic re- 
laxation rates can be used to determine the location of 
charged defects limiting the mobility in DLG. Transport 
scattering rate can be extracted from Boltzmann resis- 
tivity. The momentum relxation rate characterizes the 
broadening of single particle states and can be extracted 
from the decay of Shubnikov de Haas oscillations. A 
similar technique^ has been successfully applied to the 
analysis of the standard graphene based field effect tran- 
sistors^. Since our work establishes that the inhomoge- 
neous puddles induced by the distant charged impurities 
in Si02 can be suppressed only partially by the control 
layer screening, it is important to know where the dom- 
inant charged impurities reside and the measurement of 
the single-particle momentum relaxation (sometimes also 
called the quantum scattering) time could help locate the 
charged impurities in the DLG system. 
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Figure 1. (color online) A schematic of a typical double layer 
graphene heterostructure used in Ref. 9. Top (blue) and 
bottom (red) layers correspond to the "control" and "stud- 
ied" layers respectively with typical values of d = 4nm and 
do = 20nm considered in the text. 



Very high values of the mobility of the "studied" 
layer encapsulated in hBN that is tunable in the range 
50 — 120 x 10 3 cm 2 V~ 1 s~ 1 with the screening density 
in "control" layer n c indicate the dominant role of the 
charged defects in the elastic scattering in this system. 
Furthermore, the very high value of the observed^ max- 
imum of the gate tunable mobility suggests a weak ef- 
fect of neutral (immune to screening) elastic scattering 
mechanisms that break locally the lattice symmetry of 
graphene. Therefore these DLGs may demonstrate quan- 
tum interference behavior distinct from that observed in 
typical graphene flakes^— where strong lattice symme- 
try breaking disorder determines predominantly weak- 
localization behavior. This is especially so at relatively 
low densities n s at which the effect of trigonal warping 
of the Fermi line^S can be neglected. 

Two identical carbon atoms in the unit cell of graphene 
give rise to two degenerate flavors (valleys) of chiral Dirac 
charge carriers in its electronic spectrum confined to the 
vicinities of the two inequivalent conrens of the hexagonal 
Brillouin zone. Quantum interference of two independent 
flavors of Dirac quasiparticles is associated with weak 
anti-localization behavio r 29 ' 30 . The latter is however pro- 
tected only by the symmetry between the two sublattices 
in graphene, which can be easily broken by impurities. 
Therefore quantum interference in graphene is very sensi- 
tive to the presence of elastic scattering mechansims that 
break lattice symmetry. Generic time inversion symmet- 
ric disorder in graphene can be categorized into three 
types according to their effect on quantum interference: 
(I) potential of the charged defects that does not break 
lattice symmetry; (II) smooth random vector potential 
disorder that breaks the symmetry between the two car- 
bon atoms in the unit cell of graphene (intravalley disor- 
der); (III) atomic scale disorder that breaks lattice sym- 
metery and causes mixing of the two valleys (interval- 
ley disorder) . Random vector potential disorder leads to 
suppression of the weak-antilocalization effect whereas 
the intervalley scattering restores the weak localization 
behavior typical for two dimensional electron gas. In the 
case of high quality DLG heterostructures with very weak 
lattice symmetry breaking scattering weak antilocaliza- 
tion effect of Dirac electrons plays an important role. 

Additional screening in DLG affects only the scatter- 
ing rate due to charged defects in graphene and therefore 



changes the relative strength of intra and inter valley de- 
fects with respect to the strength of potential scatterers. 
This results in the dependence of quantum correction on 
the electron density in the "control" layer n c . In partic- 
ular, an interesting possibility arises that the intervalley 
elastic scattering, while being unimportant in determin- 
ing the semiclassical Boltzmann resistivity, becomes im- 
portant in determining the tuning of the quantum inter- 
ference correction from anti-localization to localization. 

In a wide temperature range quantum interference is 
limited by decoherence rate due to inelastic electron- 
electron collisions. The latter is strongly affected by the 
screening effect in DLG for two reasons: due to suppres- 
sion of the amplitude of the screened Coulomb interac- 
tion between charge carriers, and due to enhancement of 
the diffusion coefficient that reduces the time electrons 
spend in proximity of each other where interaction is 
the strongest. The combination of these two effects re- 
sults in strong variation of the decoherence rate with the 
screening density. This gives rise to a peculiar gate tun- 
able quantum correction to the resistivity of graphene 
in DLG devices which may change sign depending on 
the screening denisty in the "control" layer especially at 
higher temperatures. It should be emphasized that this 
gate-tunable quantum correction to the Boltzmann re- 
sistivity is manifest only when the phase breaking length 
is much larger than the elastic mean free path, i.e. at 
lower temperatures, so that the quantum interference is 
operational. Since our theory treats the quantum correc- 
tion perturbatively, it applies only in the metallic regime 
at densities well above the crossover to strong Anderson 
localization. 

The most striking consequence of the gate tunable 
quantum correction can be observed in the intermedi- 
ate density regime in which the transport in the "stud- 
ied" layer is metallic, kpl ^> 1, whence the Fermi en- 
ergy is low enough so that the trigonal warping effect 
is negligible. The quantum correction in this case is 
strongly dependent on the electron density in the "con- 
trol" layer. At sufficiently strong dephasing rate (suffi- 
ciently high temperature) this results in a gate tunable 
crossover from weak-antilocalizaiton type magnetoresis- 
tance to weak-localization type. A similar crossover that 
was temperature driven was demonstrated in graphene 
previousl y 23 ' 27 , our work demonstrates a similar tunabil- 
ity with the gate voltage in DLG systems. 

Quantum part of the resistivity also includes Altshuler- 
Aronov correction to the conductivity due to the 
electron-electron interaction ;] 31 ' 32 which typically has the 
sign enhancing the localizaiton effect and is expected to 
be suppressed by the additional screening effect which 
typically reduces the interaction strength. This does not 
have any effect on low field magnetoresistance which al- 
lows to study the qunatum interference (weak localiza- 
tion) correction separately. 

The paper is organized as follows. In section II we an- 
alyze the effect of "control" layer screening on Coulomb 
potentials of the charged defects. Sections III and IV 
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discuss classical and quantum parts of the resistivity re- 
spectively. We provide a discussion of the relevance of 
our work to experiments, particularly the data of Ref. 0| 
where the DLG system was studied experimentally in 
Section V. We conclude in Section VI. 



II. SCREENING EFFECT OF THE "CONTROL" 
LAYER 

A. Linear screening in DLG 

We start with a linear screening model for DLG het- 
erostructure. For simplicity we assume that a density rii 
of charged defects is located in a plane a distance z from 
the "studied" layer (at z = 0), which can be both above 
z > and below z < it. The "control" layer is located 
at a distance d above the "studied" layer. 

We include both intra and inter layer Coulomb interac- 
tions between electrons in DLG described by a symmetric 
matrix Uw with diagonal, U\\ — U22 — v q , and off diago- 
nal elements U12 = U21 — v q e~ qd with v q = 2ne 2 /q being 
the 2D Fourier transform of the Coulomb potential. Here 
£,£' = 1,2 correspond to the "studied" and "control" lay- 
ers respectively. We calculate the interaction energy of 
an impurity charge Z\ and a conduction electron charge 
Z2 in the "studied" layer in the presence of the screening 
effect of electron gases in both "studied" and "control" 
layers. This can be found using linked cluster expansion 
of the thermodynamic potential of the DLG. In the re- 
sulting infinite series we keep only the term proportional 
to the product of the two charges that represents the 
linear screening of the interaction potential, 



Afi(Ri-Ra) = Z X Z 2 



d<\ v q e 



(2ir) 2 e*(q) 



iq(Rl-R 2 ) 
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where q, R^ are 2D wave vector and position vector, re- 
spectively, and Z\Z2eT q ^v q is the Fourier transform of 
the bare Coulomb interaction between the charges Z\ and 
Z2 ■ For electrons we can put Z\ = Z% = 1. In Eq. ((T|) we 
introduced a dielectric function, 
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where the matrix D«/(q) takes into account the depen- 
dence of Coulomb interaction energy on the spatial sep- 
aration between the charges in z direction, 
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Here xw — (T T Pi (q, t) pe> (— q, 0)) is the density-density 
correlator which is a matrix in the layer index, £ = 
1,2, including both intra and inter layer terms. The 
the density-density correlator is renormalized by the 
Coulomb interaction matrix which is given by an infi- 
nite perturbative series 3 ^. Within the Random Phase 



Approximation the result of the resummation of the per- 
turbation series satisfies a matrix Dyson equation 35-36 , 
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where \^ is the diagonal matrix of free particle polar- 
ization operators in the two graphene layers. Solving the 
Dyson equation we arrive at the dielectric function in 
the form, 
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where x = \z — d\ — \z\ + d. The dielectric function 
demonstrates a strong dependence on both the impurity 
location z and the interlayer distance d, and in the limit 
d — > 00 approaches the monolayer form l/e{q) = 1/(1 + 

Static screening is described by the zero frequency limit 
of the free particle polarization operato r 37 i 38 , 
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where 6(t) stands for Heaviside step function, g = 4 

stands for spin and valley degeneracy combined, k F and 
Vi stand for the Fermi wave vector and the density of 
states per spin per valley at the Fermi level in the i-th 
layer, respectively. 



Non-linear screening in the electron-hole puddle 
regime 



We estimate the effect of screening by the "con- 
trol" layer on the density inhomogeneity in the "stud- 
ied" layer using the simple self-consistent approach that 
had shown good agreement with more detailed Thomas- 
Fermi-Dirac numerical calcualtions in the case of mono- 
layer graphene ^ 39 ' 40 . This estimate is based on the ap- 
proximate relation (Ep) sw (V 2 ), where the amplitude 
of the fluctuation of the local Fermi energy is equated 
to the local fluctuation in the potential energy of elec- 
trons, which reflects the local electrochemical equilib- 
rium condition. We further assume that the local Fermi 
level is given by the free particle form Ep = hvpkp^ = 
hvpy/mis, where vf and n s stands for the Fermi velocity 
of Dirac electrons and local density respectively. We also 
assume that the screening of impurity potential can be 
taken into account within linear approximation Q. We 
then end up with a self-consistent equation for the typical 
density fluctuation n* . The solution of this self-consistent 
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Figure 2. The result of the self-consistent density calcu- 
lation in the "studied" layer as a function of the density in 
the "control" layer n c for a DLG with d — 4nm, see sec- 
tion lll Bl for details, (a) Solid line corresponds to z = — 20nm, 
rii = 10 12 cm -2 . The horizontal dashed line corresponds to 
the result of the self-consistent density calculation for iso- 
lated monolayer graphene encapsulated in hBN with the same 
parameters as the solid line, (b) Solid and dashed lines cor- 
respond to z = lnm and z = 5nm respectively with the 
impurity density m = 10 11 cm -2 . 

problem is plotted in Fig. [2] as a function of the "control" 
layer denisty n c . The discrepancy between the solid black 
line in Fig. [2(a) as n c —¥ and the root mean square den- 
sity in monolayer without the screening layer under same 
conditions, dashed line, originates from the screening by 
interband transitions in the "control" layer at neutrality. 
The above analysis neglects non-linear screening effects 
in the "control" layer and therefore the saturation in the 
realistic DLG device may be at a slightly different value. 
Nevertheless Fig.[2Jand the comparison to the monolayer 
case gives a reasonable estimate of the suppression of the 
density inhomogeneity by screening. 
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Figure 3. (color online) Transport scattering length in 
the "studied" layer. Solid lines from bottom to top 
(red, green, blue and magenta) correspond to n c = 
10 9 , 10 10 , 10 11 , 10 12 cm- 2 , with z = -20nm, m = 10 12 cm" 2 . 
Dashed line corresponds to isolated monolayer graphene with 
the same parameters. For comparison the effect of a rela- 
tively small concentration of charged defects, m = 10 n cm~ 2 , 
near "studied" layer, z — — lnm is shown as dash-dotted line. 
The dotted line corresponds to kpitr = 1 where quantum and 
classical parts of the conductivity are of the same order. 



is substantially less efficient for the impurities located in 
close vicinity of the "studied" layer and for the high con- 
centration of impurities in SiC>2 substrate, see the solid 
line in Fig. 12(b). 

Note that our estimate of the density inhomogeneity 
relies on the assumption that locally electron density re- 
sembles quasiparticles in clean graphene. This is unjus- 
tified in the case of strong disorder for which kp£ ~ 1, 
see Sections and HHIJ 



C. Location of impurities in DLG 

Effectiveness of screening by the "control" layer is 
strongly dependent on the location of impurites (z) which 
can accumulate in three possible areas: Relatively small 
densities of charged defects, ni ~ 10 10 — 10 n cm _2 , can 
be located at hBN-graphene interfaces^- and on top of 
the "control" layer which can be exposed to air in re- 
alistic DLGs. A much larger concentration of charged 
impurities, n, ~ 10 12 cm -2 , is expected to accumulate 
at the surface of SiC>2 substrate that is separated from 
the "studied" layer by d ~ 20nm of hBN, see Fig. Q] 
These defects do not play an important role in transport 
in the high density regime kp\z\ > 1 due to the exponen- 
tial factors in Eqs. (JTJ and However, they dominate 
transport and density inhomogeneity at very low densi- 
ties where kp\z\ < 1 and hence the long-range nature of 
the Coulomb disorder becomes important at low densities 
near the Dirac point. 

The additional screening effect in DLG is the strongest 
for impurities located on top of the "control" layer for 
which the inhomogeneity varies by an order of magni- 
tude, dashed line in Fig. [2(b). However, the screening 



III. ELASTIC RELAXATION RATE AND 
BOLTZMANN RESISTIVITY 

A. Transport relaxation rate due to charged 
defects 

Additional screening in DLG leads to a substantial sup- 
pression of the transport scattering rate which in Born 
approximation reads, 

i = ir n * EC 1 " cos ^) Kp'IWI 2 * ( £ p - v) • (s) 

Tf r ft 
tr p , 

where <p and q = |p — p'| = 2AiF|sin§| are the scat- 
tering angle and momentum of electrons on the Fermi 
surface in the "studied" layer, and |{p'| V^|p) | 2 = |(1 + 
cosy) | Pe -29 ' 2 ' is the matrix element of the Coulomb 
potential of a defect in graphene. Fig. [3] shows the vari- 
ation of the transport mean free path approximated by 
Eq. ([S]) within experimentally accessible densities in the 
"control" layer 10 9 cm~ 2 < n c < 10 12 cmT 2 . 
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Figure 4. (color online) (a) Momentum relaxation 

rate as a function of the density in the "studied" layer 
m — W 12 cm~ 2 ,z — —20nm. Solid lines from bottom to 
top (red, blue, green and magneta) correspond to n c = 
10 9 , 10 10 , 10 11 , 10 12 cm" 2 respectively. Dotted line corre- 
sponds to kpi — 1 where Born approximation becomes in- 
sufficient, (b) Ratio of the transport and momentum relax- 
ation times. Solid lines from bottom to top (blue, green 
and red) correspond to (ru, = 10 11 cm _2 ,z = lnm), (m = 
10 11 cm~ 2 ,z = 5nm) and (m = 10 12 cm _2 ,z = — 20nm) re- 
spectively with n c — 10 12 cm~ 2 . Dashed lines correspond to 
the same parameters with n c = 10 8 cm -2 . 



B. Momentum relaxation rate: location of 
impurities 

Different spatial configurations of impurities within 
DLG may be distinguished experimentally by analyzing 
the ratio of the momentum relaxation rate to the trans- 
port scattering rate. These ratios are substantially differ- 
ent in the cases of scattering on long range potential of re- 
mote and nearby impurities. The tranport rate, Eq. ([5]), 
is mostly determined by backscattering with large mo- 
mentum transfer, and therefore decays quickly with the 
distance z between the charged defects and the "studied" 
layer. The momentum relaxation rate, given by, 



1 

Tv 



2vr 



p' 



; P '|V|p)| 2 «5(e p -e p ,). 



(6) 



includes forward scattering with small momentum trans- 
fer which is only weakly dependent on z, see Fig. Ufa). 
The ratio of the momentum relaxation rate to the trans- 
port rate is shown in Fig. Efb) for the three different 
location of impurities considered in Sec. Ill CI Note that 
Born approximation used in Eq. ([5]) is insufficient when 
kplp ~ 1, £p = vft p shown as dotted line in Fig. EJa) 
and further corrections to the elastic relaxation rate have 
to be included. In the case of remote impurities in SiC>2 
there is a peculiar regime where kp£ p ~ 1 and the Born 
approximation is insufficient yet kpl tr ^> 1 which sug- 
gests diffusive transport. Note however that the pertur- 
bative analysis presented here is quantitatively correct 
only when kp£ p 1. 



Suppression of the relaxation rate is reflected in the 
Boltzmann resistivity p(n s ,n c ) = (^-kp£ t rj shown 

in Fig.[5][a). It is the strongest near the neutrality point 
in the "studied" layer in which case both, n s -C n c and 
n s -C — ~r so that the high concentration of defects in 
Si02 substrate dominates the transport properties. In 



the formal limit 2kp' = 2y/Tm s <C «2 = 87re 2 f c , («2 is 
the Thomas-Fermi screening wavevector in the "control" 
layer) the charged impurity potential becomes effectively 
short ranged which leads to the resistivity independent 
of the density n s , 



(7) 



where x = (\z — d\ — \z\ + d). This strong-screening 
limiting resistivity, however, is unlikely to be reached, 
since within the metallic regime and accessible densities 
kp /k2 is of the order of unity. Moreover, the value of 
the resistivity at low densities n s is affected by the doping 
inhomogeneity. The resistivity maximum in the presence 
of density inhomogeneity (electron-hole puddles) can be 
estimated 2 -^ 2 - using the effective medium theory^, which 
allows averaging of the resistivity of a system split into re- 
gions with random values of resistance characterized by 
a distribution P[p(n)]. The averaged resistance in this 
case is given by the solution of a self-consistent effective 
medium equation, 



/ 



dnP[n] P{ ;\ PEMT = 0. 
p(n) + PEMT 



(8) 



We approximate the distribution of the density fluctua- 
tions by a Gaussian with zero average and variance given 
by n rms — (V 4 ), which compares well with full Thomas- 
Fermi-Dirac approach^ in the case of disordered mono- 
layer graphene. The results are shown in Fig. [SJa) as 
dashed lines. Note however that due to the relatively 
weak doping inhomogeneity induced by the remote im- 
purities in Si02 the value of the metallicity parameters 
2 ^ kp(n rms )£tr < 4 and kp£ p ~ 1 indicate that quan- 
tum corrections could be important in this low-density 
regime. However, the long-range character of the dis- 
order potential makes quantitative analysis of quantum 
interference difficult^r— in this regime. 

At higher density n s the random inhomogeneity can be 
neglected and transport is metallic, therefore the linear 
screening model and the Boltzmann resistivity provide an 
accurate descirption of the system. At even higher densi- 
ties n s irz 2 > 1 the exponential factors in Eqs. (|U [5J lead 
to the strong suppression of the effect of remote impuri- 
ties in SiC-2. As a result charged defects in the vicinity 
of graphene layer and neutral scatterers determine the 
resistivity. Finally, at the very high density kpd > 1 the 
typical monolayer behavior of the resistivity is restored, 
see Fig. [5£a)- 
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n , 10 cm 



Figure 5. (color online) Resistivity of the "studied" layer 
p(n s , n c ). (a) Solid lines from top to bottom (red, green, blue 

", 10 11 , W 12 cm' 2 
-20nm, d = 4nm. 



and magenta) correspond to n c — 10 9 , 10 10 



respectively, with 



lO^crrT 



Dashed lines correspond to inhomogeneity induced maximum 
resistivity with the same parameters from top to bottom, (b) 
Resistivity at intermediate density n s = 10 11 cm -2 . Top (red) 



line corresponds to rii 



= — 20nm, bot- 



tom (black) and middle (blue) lines correspond to n;. 
10 11 cm -2 and z — 5nm, lnm, respectively. 



D. Intermediate density regime 



There exists a regime of densities n s where Eqs. (4,5) 
provide a quantitatively accurate description of the de- 
pendence of the elastic relaxation rate on both densi- 
ties n s and n c . Formally, the constraints defining this 
regime are: (i) the effect of the density inhomogene- 



ity on the resistivity is negligible, y (n s — n s ) <C n s , 

where O stands for spatial average of O; (ii) transport 
in the "studied" layer is metallic n s ^> -W. In addi- 

tion we assume (iii) that the denisty in the "studied" 
layer is sufficiently small, n s <C n Cl so that the resis- 
tivity changes substantially with changing n c . Here n c 
is limited by the experimentally accessible density range 
n c < 5 x 10 12 cm -2 ; and (iv) the trigonal warping of 
the Fermi line has a negligible effect on weak localiza- 
tion corrections 2 ^, tJ 1 /t _1 = 2 (-K/ihTris) 2 <C 1, with 
fi = and a « 1.42Ais the distance between the 
nearest-neighbors in the hexagonal lattice of graphene. 
The latter constraint simplifies the analysis of character- 
istics of various scattering mechanisms that determine 
the quantum transport in the system, (v) Scattering 
due to Coulomb defects remains dominant fcjrz < 1. 
Throughout the text we refer to the range of densities 
n s where all constraints (i)-(v) are satisfied as the "in- 
termediate density regime" . 

The intermediate density regime can be used for a 
quantitative analysis of both Boltzmann and weak lo- 
calization contributions to the resistivity. The effect of 
screening electron density n c in the "control" layer on 
the resistivity is illustrated in Fig. [SJb). The effect of 
n c on the weak localization correction in the DLG in the 
intermediate density regime is discussed below. 



IV. QUANTUM CORRECTION 

We generalize the pcrturbative theory of weak localiza- 
tion in graphen o 29 i 47 to the DLG heterostructure in the 
intermediate density regime defined in Sec. IIIIDI above. 



A. Perturbative theory of weak localization 

We briefly review the theory of weak localization cor- 
rection to the resistivity of graphene. Weak localization 
originates from the interference contribution to the prob- 
ability to scatter backwards. The latter is given by the 
modulus square of the sum of the quantum mechanical 
amplitudes Ai associated with each possible backscat- 
tered trajectory, P ~ |^Ai| 2 . This contains the clas- 
sical part, P c i ~ l^j| 2 > an d the quantum interfer- 
ence contribution, SP ~ ^AiA*. The latter vanishes 
for generic trajectories. However, in case A* is a time 
reversal image of Ai, which is possible for loop-shaped 
self-intersecting trajectories, the phase factors cancel ex- 
actly, 5P ~ J2' l^l 2 - Therefore the average contribution 
of this special type of trajectories to the backscattering 
probability is non-zero. SP for a given self-intersecting 
trajectory is proportional to the ratio of the width of 
the trajectory VF^dt to the 2D volume Dt that can be 
spanned by diffusion during time t. The resulting cor- 
rection to the conductivity So ~ J Xvpdt/ (Dt) ~ In t v /t 
is proportional to the log of the ratio of the travel times 
along the longest t v and the shortest r trajectories and 
is typically of the order of the conductance quantum. 
The longest trajectories t v that limit the extent of quan- 
tum interference are given by either the system size or 
by phase breaking effects. 

Two degenerate flavors (two valleys) of chiral Dirac 
fermions in graphene are characterised by an isospin de- 
gree of freedom that has a fixed projection on the di- 
rection of momentum. Therefore backscattering of chiral 
charge carriers is accompanied by the revesal of isospin 
which gives rise to an additional phase difference ttN be- 
tween the two images of a looped backscattered trajec- 
tory, where N is the winding number of the trajectory. 
The interference contribution to conductivity in this case 
acquires an additional minus sign which results in anti- 
localization 30 effect in contrast to the usual weak local- 
ization in 2D 48 . In realistic devices the vector potential 
(intravalley) disorder, characterized by a phenomenolog- 
ical rate r" 1 , breaks the symmetry between the two sub- 
lattices of graphene and therefore suppresses the anti- 
localization effect. Atomic scale disorder, characterized 
by Tj~ , scatters electrons between different valleys and 
therefore restores the localization sign of the interference 
correction. Quantitatively this effect is described by a 
perturbation expansion in kp£ 3> 1 which results in the 
conductivity correction 29 , 
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Here the first term on the left hand side describes in- 
tervalley interference sensitive only to the ratio the 
second term describes the intravalley interference which 
is also sensitive to soomth intravalley scattering rate t" 1 
due to the vector potential disorder and the total trans- 
port scattering rate ^- ~ The left hand side in the 
latter equation, given by Eq. ([5]), is obtained modelling 
charged defects by a Gaussian disorder with a finite range 
in space and zero average of the potential. 

In Fig. IUJa) we show the lines in the parameter space 
that separate the negative 5a < 0, localization-like cor- 
rection from the positive one 5a > 0. In the presence 
of strong vector potential disorder, t z ~ r tr the intraval- 
ley interference term is strongly suppressed. This results 
in the weak localization sign of the correction even in 
the case of relatively weak intervalley scattering, see the 
blue line in Fig EJa). By contrast, in the presence of 
relatively weak vector potential disorder T-T 1 <C t^, 1 , a 
substantially stronger intervalley scattering is required 
to compensate for the intravalley term in @ and for 
the correction to have the localization sign 49 . We em- 
phasize that these analytical considerations are entirely 
restricted to the perturbative regime, kpf-p 3> 1 of the 
quantum interference correction to the conductivity and 
do not apply to the strong localization regime of metal- 
insulator transitions. 



B. Intervalley scattering 

Microscopic origin of the intervalley scattering at in- 
termediate densities therefore deserves a more careful 
consideration. Electrostatic potential can give rise to a 
strong intervalley scattering only if its amplitude is larger 
than the bandwidth in which case the lowest order Born 
approximation is insufficient, see Appendix and Ref. [5CJ 
for details. Such strong potentials can be induced by lat- 
tice defects and adsorbate molecules chemically bound to 
carbon atoms of the graphene lattice. Scattering charac- 
teristrics of such defects depend strongly on the detailed 
nature of the electronic structure, type of the chemical 
bonding and the position of the defect in the unit cell 
of graphene^. Such details are beyond the scope of the 
present analysis. Instead we employ a more generic phe- 
nomenological model of two types of scatterers namley 
resonant scatterers and short range disorder. The for- 
mer describes the effect of monovalent chemical groups 
and vacancies in the graphene lattic o 52 ' 53 . The scattering 
rate in this case decreases with density 1 ~ ( s/n In n) 1 
which in monolayer graphene is practically indistinguish- 
able from that due to charged defects ~ 1/y/n, see 
Ref.0 Short range scatterer model describes atomic de- 
fects that do not give rise to resonant scattering. The 
corresponding scattering amplitude grows with density 
t _1 ~ %fn and the resulting contribution to the conduc- 
tivity is density independent. In the case of both sources 
of intervalley scattering a relatively weak r" 1 is sufficient 
for the quantum correction to have the weak localization 



sign, FigEHa). 

There exists substatnial experimental evidence for the 
presence of intervalley scattering in graphene. In parti- 
cluar, the observation of negative weak localization mag- 
netoresistance suggests a rather high intervalley scatter- 
ing rate corresponding to Lj„ = \JDt^ ] ~ 0.1 — 0.3/zm at 
high densities n > 10 11 cm~ 2 ^ ~ 28 ' 54 . Alternative mea- 
surements of a weak D-peak in the Raman spectrum 
demonstrated the presence of intervalley scattering at 
high energies ~ leV— . More detailed quantitative char- 
acteristics of the intervalley scattering such as density 
dependence are not currently available experimentally. 



C. Phase breaking effect 

The size of the self-intersecting trajectories that give 
rise to the interference correction is limited by either the 
system size or the phase breaking effect. The latter can 
be caused by scattering on magnetic defects^ or inelastic 
electron-phonon or elelctron-electron collisions^. Both 
latter effects result in temperature dependent dephasing 
rate t" 1 , which is linear t^ 1 ~ T in the case of electron- 
phonon scattering and linear or quadratic in the case of 
electron-electron interactions in the diffusive and ballis- 
tic regime respectively. Experiments report typical val- 
ues for the dephasing length ^ ' Dt v ~ 2/xm at the lowest 
temperatures T < IK—"— and L v ~ 0.2/xm at higher 
temperatures T ~ 2QK. More recent measurements on 
large area epitaxial graphene devices accessed very high 
temperature regime T ~ 1QQK corresponding to a short 
dephasing length L v ~ 0.04/i m 56 i 58 . The available exper- 
imental information about the density, disorder, and tem- 
perature dependent phase breaking length in graphene is, 
however, not extensive, and more work is needed before 
definite quantitative statements can be made in details. 



D. Effect of screening on the dephasing rate due to 
inelastic electron-electron collisions 

In this section we estimate the effect of additional 
screening on the dephasing rate caused by inelastic 
electron-electron collisions with low energy transfer hu < 
ksT. We estimate the phase breaking rate up to a numer- 
ical coefficient using the self-consistent approac h 31 ' 32 ' ^° 
which is sufficcnt to analyze the dependence of the phase 
breaking rate on the densities in the two layers of DLG. 
There are two distinct dephasing temperature regimes, 
diffusive, ksT » hr^ 1 , and ballistic, ksT -C frr^ 1 
which correspond to distinct temperature dependences 
ofr" 1 . 

In the diffusive regime, fcsTr tr -C h, electrons un- 
dergo diffusive motion which increases the time they 
spend in vicinity of each other (where interaction is the 
strongest) and therefore enhances the inelastic scattering 
rate. Electron-electron collision rate is given by Fermi 
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Figure 6. (color online) (a) Weak antilocalization to locaiza- 
tion crossover boundary in the parameter space of disorder 
in DLG. Magenta, green and blue lines correspond to the 
fixed values of t v /t z — 100, 10, 3 respectively. Red line cor- 
responds to monovalent adsorbate disorder. The dashed lines 
show example trajectories that a DLG initially at point A 
takes with changing n c in case of fixed t v , straight line, and 
coherence limited by the system size, the curved trajectory 
~ \J Ttpjrtr. Dash dotted line corresponds to the trajectory 
with coherence limited by electron-electron interactions in the 
diffusive regime Eq. (|13|). (b) Weak localization correction 
for different values 0.1/^m < upn < lfim from bottom to 
top with n B = 5 x 10 10 cm -2 . We assume rii — I0 12 crn~ 2 at 
z = — 20nm and dephasing is limited by the electron-electron 
interactions in the diffusive regime Tr <C 1, Eq (113[) . We took 
vft v = 0.3/im at n c = 10 10 cm -2 . 



golden rule, 



1 



— = Akv 6 



1 dujde'F{e,e',uj)W 2 (Lu), (10) 

to 

F(e,e',Lj) = f' e (l ~ / e _ u )(l - fc+u) + (1 - ftfe-vU+v, 



where f e = 1/(1 + exp(^e)) and /T 1 = k B T. 
action kernel in (|10[) is given by, 



The inter- 



(ii) 
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es). 



The matrix element of the screened interaction is written 
in terms of exact eigenstates \a) of non-interacting Dirac 
Hamiltonian in the presence of a given disorder poten- 
tial. The screened electron-electron interaction potential 

which includes the 



in DLG is given by U u (q, w) = j^jy 
dynamical dielectric function of Eq. Q where we need 
to put z = for interaction within the "studied" layer. 

The expression for the collision rate (fT0|) has to be 
averaged over disorder realizations. The result of this 
averaging can be then expressed in terms of the diffu- 
sion propagator in graphene D{q 1 uj) which describes the 
transport of the charge density through the disordered 
system. The diffusion propagator is found by solving the 
diffusion equation (<9 t — ADi) D(r, r') = S(r — r') in the 
"studied" layer—. Here Di is the diffusion coefficient in 
the "studied" (i — 1) or "control" (i = 2) layer respec- 
tively. As a result the Eq. (fTU)) at ql\ <C 1 and u)T tr <C 1 



takes the form, 



1 



8// 



cLucPq w\U(q, lo)\ 2 



(2nf 



sinht 



[nD(q,uj)} 2 (12) 



where we keep only the logarithmically divergent terms. 



Also we introduce a low energy cutoff ui 



which 



reflects the fact that only collisions with large enough 
energy transfer uj > r^ 1 produce a substantial dephasing 
of the electrons with diffusion time limited by T yj 31 i 32 i 62 . 
Solving (fT2"|) self-consistently gives an estimate for the 
dephasing rate. 

At high densities such that ^ > 1, i — 1,1, the 
screening is ballistic for which case we can approximate 
qv q Xu(q,0) ~ K t and Eq. flT2J| gives, 



T 



with 



- i c 2 

2 ' 1 



(13) 



(14) 



where we introduce the screening parameter Ki = 87re 2 i^ 
for i-th layer and £ = 1 + 2dK 2 - 

At lower densities the polarization operator aqcuires a 
diffusion pole, 



Xifew) « gvi 1 + — 



1 



— iu) + Diq 2 



(15) 



The presence of the diffusion pole affects only the coef- 
ficient in Eq. (fT5|) . In the case of short screening length 
Ki = %-Ke?Vi > we obtain for the coefficient C — >• C, 



C 



c + f ) + 



K 2 ( ] 



«2 Eh 



tt(C+f) (l+(C + 2f) £) 



(16) 



In the limit of K2/K1 3> 1 the results (jT6j) and (jT4j) both 
approach, 



G 



2 k| 



(17) 



The dephasing rate p^|) is strongly dependent on the 
screening density n c in the "control" layer through the 
diffusion coefficient D\ and the factor C or C. 

In the high temperature Tr fr 3> 1 ballistic regime the 
polarization operator is approximated by the free particle 
fo rm 3L3£, Xl ( g , w ) w + a t w/(u F ff) < 1. 

The dephasing rate in this case is given by^ 3 -, 



(fc B T) 2 



16^i 



c + — 



hi 



k B T 



(18) 



where parameters Ki = S"ne 2 Vi and Q = 1 + 2d,K2 have the 
same meaning as above and we assume Kid <C 1. 
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E. Gate tunable quantum interference correction 



In the presence of mixed scattering sources the sup- 
pression of the transport scattering rate with growing 
n c results in the increased relative strength of neutral 
scattering mechanisms. This is a simple consequence of 
the Coulomb scattering being effectively screened out by 
the "control" layer. As a result the amplitude of the 
quantum interference correction varies as a function of 
n c . To illustrate this we assume a large concentration of 
charged defects rn located in the Si02 substrate and a 
small concentration of resonant scatterers Uj. Then, we 
have Tj" 1 s» t^ 1 + 2r~ 1 where we use t~ x s» r^" 1 valid 
for a small concentration of single site defects. In the 
present analysis we ignore the effect of screening by the 
"control" layer on the shift of the chemical potential in 
graphene due to charged defects. 

Depending on the mechanism limiting the phase co- 
herence a variety of behaviors of the conductivity cor- 
rection as a function of the screening density n c is pos- 
sible, see Fig. EJa). For screening independent r^ 1 due 
to dephasing by magnetic impurities or electron-phonon 
interactions the system initially at point " A" in Fig. Eta) 
takes a trajectory represented by the straight dashed 
line. In smaller devices coherence is more likely to be 
limited by the system size, in which case the effective 
rate t^ 1 = = "|^ r , depends on the diffusion coeffi- 
cient and hence is sensitive to the screening effect of the 
"control" layer. This regime corresponds to the curved 
dashed trajectory in Fig. Eta). The effect of screening 
by the "control" layer is the strongest in the case of de- 
phasing limited by electron-electron interactions. In this 
case the ratio r^/rtr depends on n c through the coef- 
ficient C in Eq. (TT51) . and the system takes the trajec- 
tory shown as the dash-dotted line in Fig. |6ja) . This is 
the situation most likely to be realized in experiments 
on large exfoliated flakes of graphene. Dependence of 
the weak localization correction on n c for this case is 
given by Eq. (J9j) combined with ([5]) and (|T6l) . which is 
shown in Fig. E^b). The effect of additional screening by 
"control" layer therefore expands the parameter range in 
which weak-localization sign of the quantum correction 
could be observed, and thus the possibility of a transi- 
tion from the weak anti-localization to weak localization 
as a function of the "control" layer density (i.e. gate- 
tunable) in the nominally metallic studied layer is very 
high in this case. This predicted transition has not yet 
been seen in experiments, but may actually exist in the 
samples of Ref. |9| in the metallic regime of the studied 
layer. 

Weak field magnetoresistnace allows an unambigous 
measurement of the weak-localization correction. At 
small fields B <C B v — h/e(Dr ip ) the magneto-resistance 



can be approximated by, 




In this regime we expect a gate tunable crossover from 
anti-localization to weak localization magnetoresistance 
to be observable especially at higher temperatures. 

V. COMPARISON WITH EXPERIMENT 

In this section we discuss the relation of the theory 
presented above to the recent transport measurements 
in DLG—. The clever DLG design of Ref. |9j provides an 
opportunity to use the additional gate tunability of the 
resistivity in the "studied" layer, 

p(n s ,n c ) = p c (n s ,n c ) + p R (n c ), (20) 

to distinguish the contributions of charged pc{n Sl n c ) and 
neutral Pr{ti c ) scattering mechanisms in the intermedi- 
ate denisty regime, see Sec. IIIID1 The latter can be de- 
duced from the dependence of the mobility of the "stud- 
ied" , yu(n c ), on the "control" layer density reported in the 
Supplementary Material of Ref. 0. On the one hand these 
data are consistent with the presence of a high concen- 
tration of charged defects rii ~ 10 12 cm~ 2 at the surface 
of SiC>2 substrate located at z = — 20nm away from the 
"studied" layer-. On the other hand a similar behavior 
could result from a smaller concentration of charged de- 
fects rii « 10 11 crn -2 in close vicinity of the "studied" 
layer or on top of the "control" layer, see Fig. [5jb) for 
comparison. The two situtations may be distinguished 
experimentally by comparing the transport and momen- 
tum relaxation rates, see Fig. EJb). At low carrier den- 
sities however the precise location of impurities does not 
play an important role and therefore the substantially 
larger concentration of defects likely to be present in Si02 
dominates the transport properties. 

The situation is different in the case of density inho- 
mogeneity. We find that DLG is relatively inefficient in 
screening out electron-hole puddles induced by the high 
concentration of defects in Si02 as well as defects in close 
proximity to the "studied" graphene layer itself. Both 
allow some "control" in the expreiments, the former can 
be reduced by changing the substrate and the latter may 
be reduced by improving the sample preparation. Mor- 
ever, both improvements were recently implemented in 
suspended devices^. 

In the absence of more systematic data we can estimate 
the residual mobility associated with neutral scattering 
mechanisms from the maximum of the mobility curve, 
p(n c ) at the maximum density in "control" layer reported 
in the Supplementary Material of Ref. 0. The result 
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is fi(n s — 10 11 ) ~ 1.2 x 10 5 cm~ 2 V~ 1 s~ 2 , which corre- 
sponds to a relatively small density of resonant scatter- 
ers, such that vfTr > 0.8fim suggesting that the DLGs 
analyzed in Ref. [9| may demonstrate gate tunable weak 
anti-localization to localization crossover at intermediate 
densities. 

Finally, we discuss the effect of coupling between 
graphene and hBN substrate which gives rise to spec- 
tacular new features in the electron spectrum at high 
densities^. However, at low carrier densitites this ef- 
fect is not expected to be dramatic. Transport charac- 
teristics of graphene could be affected in three possible 
ways: (i) triangular symmetry of hBN breaks the sub- 
lattice symmetry of graphene and therefore gives rise to 
an additional trigonal warping of the Fermi surface^; 
our estimate shows negligible effect on the weak localiza- 
tion correction at low densities n s < 10 11 cm~ 2 ; (ii) hBN 
induced Moire pattern may dominate the density fluc- 
tuation at the neutrality point. This effect however has 
not been clearly identified in the local probe data^ to this 
point and therefore we expect it to be smaller than cur- 
rent levels of inhomogeneity induced by charged defects, 
(iii) Additional intervalley scattering may be caused by 
strong coupling of graphene to hBN. However, existing 
experimental data and ab initio calculations suggest a 
rather weak coupling which is not expected to give rise 
to strong intervalley sattering. Nevertheless, a more de- 
tailed analysis is required to completely rule out all of 
these possibilities. 



VI. DISCUSSION AND CONCLUSION 

Our goal in this work has been a careful theoret- 
ical study of the screening by the "control" layer in 
determining transport properties of the studied layer 
in double-layer graphene systems experimentally intro- 
duced in Ref. U The problem of interest to us is both 
subtle and complex, involving many disparate aspects of 
electronic transport phenomena, and our study is at best 
an approximate one because of the highly complex na- 
ture of the problem. In particular, we are interested in 
studying both semiclassical Boltzmann and quantum lo- 
calization contributions to the DLG resistivity, including 
the "control" layer screening effect, on an equal footing as 
much as possible. Since potential, intravalley and inter- 
valley elastic scattering affect quantum localization prop- 
erties qualitatively differently, our work must include all 
mechanisms in the presence of "control" layer screening 
with equivalent considerations for all scattering processes 
contributing to the semiclassical and the quantum part 
of graphene resistivity. Since the phase breaking length 
arising from inelastic scattering processes is an impor- 
tant ingredient in determining localization effects arising 
from quantum interference, we consider electron-electron 
interaction induced inelastic phase decoherence effects in 
our theory, neglecting electron-phonon interaction which 
is known to be weak in graphene^. We also neglect all 



direct contributions to transport from inelastic scattering 
processes since our interest is relatively low-temperature 
transport where quantum interference may play a role 
rather than high-temperature transport where inelastic 
phonon scattering plays a role. We investigate the role 
of the "control" layer screening on the density inhomo- 
geneity (i.e. electron-hole puddles) near the Dirac point 
finding that screening by the "control" layer is ineffec- 
tive in suppressing the electron-hole puddles arising from 
charged impurities in the Si02 substrate far from the 
"studied" graphene layer (in fact, this is an important 
qualitative new result of our work). Finally, our work 
treats both long-range Coulomb disorder arising from 
random charged impurities in the environment and the 
short-range disorder arising from neutral atomic defects 
on an equal footing. 

Since the theory involves (at least) six different inde- 
pendent characteristic length scales (elastic mean free 
path due to long-range and short-range disorder and 
intra- valley and inter-valley scattering length, the in- 
elastic phase breaking length, the Fermi wavelength, 
the density fluctuation correlation length, the thermal 
length, etc.), it is necessarily a complex problem neces- 
sitating various approximations focusing on different as- 
pects of the transport phenomena. We concentrate on 
low-temperature transport at relatively high "metallic" 
densities where the quantum interference contribution to 
the resistivity can be treated as a perturbative weak lo- 
calization correction to the semiclassical Boltzmann re- 
sistivity. This precludes us from commenting directly on 
the very low density nonperturbative insulating behav- 
ior (as observed in Ref. 0) from the perspective of strong 
Anderson localization phenomenon, but our theory es- 
tablishes the clear possibility of a gate-tunable transition 
from weak anti-localization to weak localization behavior 
in graphene at metallic densities, which would be a neces- 
sary precursor to a possible strong localization crossover 
behavior at low densities. Whether such a weak local- 
ization transition, which must precede any strong local- 
ization induced low density metal-insulator transition, is 
operational in Ref. IH was unfortunately not studied there, 
and the prediction of such a gate-tunable transition in the 
weak localization behavior is an important new predic- 
tion of our work. We emphasize that such transition must 
be studied at relatively low temperatures, and it may 
be missed at the higher temperature range (10-100K) of 
study used in Ref. 0- In addition, the strong localiza- 
tion induced insulating behavior at low carrier density 
should be exponential in temperature, and not a power 
law as observed in Ref. |g. Given our concrete theoretical 
predictions, we hope that future experiments on DLG 
systems will resolve the important question of a possi- 
ble low-density quantum localization transition in clean 
graphene. Our work shows that it is possible for the in- 
tervalley scattering in graphene to be strong enough to 
induce localization, but at the same time be weak enough 
not to strongly suppress high-density mobility, particu- 
larly in the presence of screening by the "control" layer. 
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In conclusion, we analyze the effect of gate tunable 
screening provided by the additional "control" layer in 
DLG on both classical and quantum parts of the resis- 
tivity of a high quality "studied" layer in which elastic 
scattering is limited by charged disorder. We find that 
(i) the additional screening in DLG is relatively ineffi- 
cient at screening out electron-hole puddles induced by 
the high concentration of defects in SiC>2 as well as de- 
fects in close proximity to the "studied" graphene layer 
itself. In the latter cases, the screening provides roughly 
a factor of 2 suppression within the "control" layer den- 
sity range of 10 10 — 10 12 cto~ 2 . Therefore the suppression 
of electron-hole puddles itself is insufficient to explain 
the contrast between the strongly insulating behavior ob- 
served in DLG 9 structures and the absence of such in 
stand-alone graphene samples including highly homoge- 
neous suspended samples^, (ii) The screening effect of 
"control" layer results in a relatively strong suppression 
of the elastic relaxation rate resulting in gate tunable 
mobility of the "studied" layer in the DLG at interme- 
diate density, n s . (iii) We show that the location of the 
charged defects limiting the modbility of the DLG can be 
determined by measuring the ratio of the momentum re- 
laxation and transport relaxation rates, (iv) Additional 
screening effect of the "control" layer strongly suppresses 
the potential of Coulomb interaction between charge car- 
riers, (v) Combined effect of (ii) and (iv) results in a 
strong suppression of the dephasing rate due to inelastic 
electron-electron collisions which determines the sign and 
the magnitude of the quantum interference correction to 
the resistivity, (vi) As a consequence of (ii) and (v) weak 
localization correction to the resistivity is gate tunable. 

(vii) The very low concentration of resonant or short 
range scatterers that has little effect on the Boltzmann 
part of the resistivity in high quality "studied" layer in 
DLG nevertheless can provide sufficient intervalley scat- 
tering for the quantum correction to resistivity to have 
a negative (weak localization) sign. Moreover, the addi- 
tional screening in DLG improves the coherence in the 
system and therefore expands the range in which quan- 
tum correction of weak localization sign is observable. 

(viii) We find that high quality DLG structures where 
charged defects are the dominant source of elastic scat- 
tering manifest a peculiar gate-tunable crossover between 
weak anti-localization and weak localization magneto- 
resistance. The latter could be observed in the inter- 
mediate density regime in which electron density in the 
"studied" layer is sufficiently high for random inhomo- 
geneity to be neglected and the transport to be metal- 
lic and at the same time is low enough for the effect of 
the "control" layer screening to be substatntial, see also 

secinrDj 

In this paper we developed a theoretical framework 
that can be used for a systematic experimental charac- 
terization of the symmetry breaking mechanisms in DLG 
in the intermediate density regime which could allow 
further understanding of the observed strongly insulat- 
ing behavior at low densities. This work could also be 
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Figure 7. The ratio of intravalley scattering time (n) to in- 
tervalley scattering time (T2), ri/ra, or equivalently, the ratio 
of the scattering rates, 72/71, where ji — h/n for the long 
range Coulomb potential, (b) and (c) show n and t%, re- 
spectively, as a function of density for an impurity density 



10 11 cm 2 . Note that n oc n 1 ^ 2 and ti 



-1/2 



used to describe a more generic double layer heterostruc- 
ture involving materials other than graphene that realize 
other universality classes of Dirac electrons such as sur- 
faces of strong and weak topological insulators. Such 
heterostructures could be useful to study experimentally 
gate tunable quantum interference and Anderson local- 
ization effects in these systems which are a topic of active 
current theoretical interest^. 
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Appendix: Intervalley scattering due to charged 
defects 



Intervalley scattering rate due to Coulomb impurities 
can be separated from the intravalley part by taking low 
and high momentum parts of the sum over p' = p + q + 
G in Eq. ([5]). Where G is the reciprocal lattice vector. 
For intravalley scattering we include only G = which 
results in, 



l mv 
IT 



r? 



2k r 



q 2 dq 

p 
F 



V(q) 



y/1 - {q/2k F f (A.l) 



For intervalley scattering we take the smallest allowed 
reciprocal vector that connects the two corner of hexag- 
onal Brillouin zone, |G| = |K — K'| = 2ir/3a, where 
a = 1.42Ais the C-C distance. Then the intervalley scat- 
tering time is given by 



1 TliV 

7c ~ 



1 



h fc|G-k| 



1+ 



qdq[V(q)} 2 sin ( 



(A.2) 



where q± = | G — p| ±p and sin 6 is given by the following 
relation 



- jr 



2|G - p|pcos(9. 



(A.3) 
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In Fig. [7] we show the calculated ratio of intravalley scat- 
tering time to intervalley scattering time, r-p /t§ as a 
function of density for the long range Coulomb impuri- 
ties. Since the ratio is very small the intervalley scatter- 



ing plays a little role in the transport when the Coulomb 
disorder dominates. At the same time a phoneomenolog- 
ical model of delta-range disorder gives ti/t 2 = 1 at all 
densities. 
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